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Abstract 

Background: lA^/?/(/ transcription factors are one of tine largest families of transcriptional regulators in plants. WRKY 
genes are not only found to play significant roles in biotic and abiotic stress response, but also regulate growth 
and development. Grapevine {Vitis vinifera) production is largely limited by stressful climate conditions such as cold 
stress and the role of WRKY genes in the survival of grapevine under these conditions remains unknown. 

Results: We identified a total of 59 VvWRKYs from the V. vinifera genome, belonging to four subgroups according 
to conserved WRKY domains and zinc-finger structure. The majority of VvWRKYs were expressed in more than one 
tissue among the 7 tissues examined which included young leaves, mature leaves, tendril, stem apex, root, young 
fruits and ripe fruits. Publicly available microarray data suggested that a subset of VvWRKYs was activated in response to 
diverse stresses. Quantitative real-time PGR (qRT-PGR) results demonstrated that the expression levels of 36 VvWRKYs 
are changed following cold exposure. Gomparative analysis was performed on data from publicly available microarray 
experiments, previous global transcriptome analysis studies, and qRT-PGR. We identified 15 VvWRKYs in at least two of 
these databases which may relate to cold stress. Among them, the transcription of three genes can be induced by 
exogenous ABA application, suggesting that they can be involved in an ABA-dependent signaling pathway in response 
to cold stress. 

Conclusions: We identified 59 VvWRKYs from the V vinifera genome and 15 of them showed cold stress-induced 
expression patterns. These genes represented candidate genes for future functional analysis of VvWRKYs involved in the 
low temperature-related signal pathways in grape. 
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Background 

Plants have a variety of defense mechanisms to protect 
themselves from adverse environmental effects. Families 
of transcription factors are involved in these processes 
by functioning to reorganize gene expression patterns. 
The WRKY family is among them and plays key roles in 
modulating genes expression during plant defense in re- 
sponse to pathogens [1,2]. The WRKY transcription fac- 
tors were first identified in sweet potato (SPFl) as DNA 
binding proteins [3]. Two similar genes {ABFl and 
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were found in wheat during germination [4]. Subse- 
quently, Rushton et al. [5] reported the identification and 
characterization of WRKYh WRKY2 and WRKY3 from 
parsley {Petroselinum crispum) and proposed these genes 
belong to a gene family. This gene family was named 
WRKY due to a conserved region (WRKYGQK) that was 
identified in the N-terminal amino acid sequence of all the 
members [4,5]. Further studies showed that the conserved 
WRKY domain had other forms such as WRKYGKK and 
WRKYGEK [6], or the WRKY domain could be replaced 
by WKKY, WKRY, WSKY, WIKY, WRIC, WRMC, WRRY 
or WVKY [7,8]. 

According to variation in WRKY domain and a zinc 
finger motif in the C-terminus, WRKY proteins were di- 
vided into four groups [9,10]. WRKY proteins with two 
WRKY domains composed group I. Groups II and III 
were characterized by a single WRKY domain. Group II 
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WRKY proteins were further subdivided into five or 
more subgroups based on short conserved structural 
motifs while group III proteins contained a variant zinc- 
finger which ends with HXC. Finally, group IV WRKY 
proteins contained the WRKY domain, but lack a com- 
plete zinc-finger structure in the C-terminus. WRKY 
proteins usually functioned as transcriptional regulators 
via binding to W-boxes (TTGACC/T) in the promoter 
regions of down-stream genes and clusters of W-boxes 
had an amplified effect [3-5,11-15]. However, some other 
studies have found that some WRKY proteins bind to 
the PRE4 element (TGCGCTT), SURE element (TAAA 
GATTACTAATAGGAA) or SURE-like element and the 
WK box (TTTTCCAC) [2]. 

WRKY proteins have been found to play essential roles 
in pathogen defense in response to bacteria [16,17], 
fungi [18,19], and viruses [20,21]. Evidence also sup- 
ported that WRKY transcription factors were involved in 
modulating gene expression in plants during abiotic 
stresses such as cold [22,23], salt [24,25] and drought 
[26-28]. Besides roles in response to biotic and abiotic 
stress, WRKY proteins were also implicated in processes 
that modulate plant developmental processes such as 
morphogenesis of trichomes and embryos, senescence, 
dormancy, and metabolic pathways [2]. 

Grape is one of the most important fruit crops world- 
wide. The productivity of grapevines is largely limited by 
disease pressure and stressful fluctuations in environ- 
mental conditions. Due to their essential role in the early 
response to pathogens and abiotic stresses, several WRKY 
genes were intensively studied in grape. VvWRKYl and 
VvWRKY2, isolated from grape {V. vinifera cv. Cabernet 
Sauvignon) berries, were found to potentially participate 
in defending against fungal pathogens [18,29]. VvWRKYl 
was found involving in enhanced protection against Botry- 
tis cinerea by transactivating the VvLTPl promoter [30], 
and VvWRKY2 may regulate lignification and response to 
biotic or abiotic stresses in grapevine [31]. VpWRKYl and 
VpWRKY2, isolated from Chinese wild V, pseudoreticu- 
lata, may contribute to resistance to powdery mildew 
(Erysiphe necator) and tolerance to salt and cold stresses 
in grape [32]. VpWRKYS was found to be involved in 
pathogen defense and also interact with the salicylic acid, 
ethylene, and abscisic acid signal pathways [33]. Transge- 
netic Arabidopsis plants expressing VvWRKYl 1, isolated 
from 'Beifeng', an interspecific cultivar of V, thunbergii x 
V, vinifera, showed increased dehydration tolerance [34]. 
Its homologous gene, VpWRKYl 1, was found to serve as a 
negative regulator of disease resistance [35]. Although sev- 
eral individual WRKYg^nes have been identified in grape- 
vine, the WRKY gene family in grapevine remains wholly 
uncharacterized. 

Based on our previous transriptome analysis, we found 
that some WRKY genes respond to cold stress in different 



patterns in V amurensis (a cold hardy grapevine species) 
and V vinifera cv. Muscat Hamburg [36]. VvWRKYM 
(GSVIVT01015952001) and VvWRKY12 (GSVIVTOlOl 
2682001) were found up-regulated over 30 fold in V, 
amurensis after being subjected to cold stress but up- 
regulated to a lesser extent in V vinifera. In contrast, the 
expression of VvWRKY43 (GSVIVT01030258001) was up- 
regulated in V, vinifera (26 fold) while expression 
remained low in V. amurensis. These different gene ex- 
pression patterns in response to cold stress may be con- 
tribute to the distinctive cold hardiness between the two 
species. To further characterize how WRKY genes respond 
to freezing stress of grapevine, we initiated this study to 
identify the entire WRKY gene family in grapevine based 
on the published 12x V, vinifera cv. Pinot noir (PN40024) 
genome sequences [37]. A phylogenetic tree was con- 
structed for identified WRKY proteins and the gene ex- 
pression patterns in different tissues of V, vinifera were 
detected by RT-PCR. WRKY genes responding to biotic 
and abiotic stresses were cross-evaluated by using public 
gene-chip databases. Additionally, real time RT-PCR was 
used to detect the expression level of VvWRKYs under 
cold treatment and exogenous ABA. A comparative ana- 
lysis was conducted to identif)^ VvWRKYs that may partici- 
pate in cold signal transduction pathways in V vinifera 
using microarray data in public databases, our previously 
reported transcriptome data and qRT-PCR analysis con- 
ducted in this study. 

Results 

Identifying of WRKY transcription factors in V. vinifera 
genome 

A total of 64 transcripts in the V, vinifera genome se- 
quence were identified as possible members of the 
WRKY family. Five transcripts were excluded due to a 
lack of the conserved WRKY domain in the predicted 
amino acid sequences. The remaining 59 transcripts 
were named from VvWRKYl to VvWRKYS9 according 
to their order in the V, vinifera genomic sequence 
(Table 1). As for the previously published six WRKY 
proteins in grapes [18,29-35], each amino acid sequence 
was downloaded and BLASTp was used to find its corre- 
sponding WRKY loci in the V. vinifera genome. 

The putative genome location of each VvWRKY in the 
grape genome was shown in Additional file 1: Figure SI. 
Fifty-eight of the VvWRKYs could be mapped to 18 of 
the 19 grape chromosomes, with no VvWRKYs found on 
chromosome 3. VvWRKY4 was putatively located on the 
'Chromosome Unknown'. WRKY transcription factors 
were not evenly distributed across the chromosomes of 
the grape genome. There were most abundant on 
Chromosome 4 (8 VvWRKYs) and chromosome 7 (7 
VvWRKYs) and least abundance on Chromosome 5 and 
18 (1 VvWRKY), 
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Table 1 Identified WRKY genes in 12x V, vinifera 'Pinot Noir' genome 



Gene ID 


Gene symbol 


Subgroup 


Chromosome no. 


Peptide length 


GSVI\^01 000752001 


VvWRKYOl 


lid 


ehr7 


285 


GSVI\^01 001 286001 


VvWRKY02 


IV 


chr2 


106 


GSVI\^01 001 332001 


VvWRKY03 




ehrl_randonn 


436 


GSVI\^0 1007006001 


VvWRKY04 


1 


ehrUn 


551 


GSVI\^0 1008046001 


VvWRKYOS 


lib 


ehrl7 


606 


GSVI\^01 008553001 


VvWRKYOe 


lie 


ehrl7 


152 


GSVI\^01 009441 001 


VvWRKY07 


lid 


ehrl8 


320 


GSVI\^01 01 0525001 


VvWRKYOS 


lie 


ehrl 


190 


GSVI\^01 01 1356001 


VvWRKY09 


lib 


ehrl4 


503 


GSVI\^01 01 1472001 


VvWRKYlO 




ehrl4 


890 


GSVI\^01012196001 


VvWRKYl 1 


lie 


chrl 


284 


GSVI\^01 01 2682001 


VvWRKY12 


lib 


ehrlO 


511 


GSVI\^01 014854001 


VvWRKYlS 




ehrl9 


623 


GSVI\^01 01 5952001 


VvWRKYM 


lla 


ehr9 


279 


GSVI\^01 01 8300001 


VvWRKYlS 


lie 


ehrl5 


229 


GSVI\^01019109001 


VvWRKYie 




chr4 


487 


GSVI\^01019419001 


VvWRKYl 7 


lie 


chr2 


324 


GSVI\^01019511001 


VvWRKYlS 


III 


chr2 


343 


GSVI\^0 1020060001 


VvWRKYl 9 


lib 


chrl 


595 


GSVI\^0 1020864001 


VvWRKYlO 


lie 


ehrl2 


312 


GSVI\^01 021 252001 


VvWRKYll 


lie 


ehrlO 


279 


GSVI\^01 021 397001 


VvWRKY22 


lie 


ehrlO 


320 


GSVI\^01 021 765001 


VvWRKY23 


lie 


ehrlO 


422 


GSVI\^0 1022067001 


VvWRKY24 


lid 


chr7 


281 


GSVI\^01 022245001 


VvWRKY25 


lie 


chr7 


194 


GSVI\^01 022259001 


VvWRKY26 


lie 


ehr7 


227 


GSVI\^0 1023600001 


VvWRKY27 




ehrll 


500 


GSVI\^0 1024624001 


VvWRKY2S 




ehr6 


571 


GSVI\^01 025491 001 


VvWRKY29 


IV 


chr6 


122 


GSVI\^01 025562001 


VvWRKY30 




ehr8 


439 


GSVI\^01 026965001 


VvWRKYSl 


lie 


ehrl5 


349 


GSVI\^0 1026969001 


VvWRKY32 


lie 


ehrl5 


202 


GSVI\^0 1027069001 


VvWRKY33 


III 


ehrl5 


361 


GSVI\^01 0281 29001 


VvWRKY34 


lie 


chr7 


243 


GSVI\^01 0281 47001 


VvWRKY35 


lie 


chr7 


303 


GSVI\^0 1028244001 


VvWRKY36 


lib 


ehr7 


480 


UoVIV 1 U 1 Uzo/ 1 oUU 1 


VVVVnrK Yj/ 


III 
ill 


enr i o 


jOj 


GSVI\^0 1028823001 


VvWRKY3S 


lie 


ehrl6 


183 


GSVI\^01 029265001 


VvWRKY39 


lid 


ehrll 


280 


GSVI\^0 1029688001 


VvWRKY40 


lib 


ehrl2 


491 


GSVI\^0 1030046001 


VvWRKY41 




ehrl2 


365 


GSVI\^01 0301 74001 


VvWRKY42 


III 


chrB 


332 


GSVI\^01 030258001 


VvWRKY43 




ehr8 


514 


GSVI\^01 030453001 


VvWRKY44 


lib 


ehrl2 


499 



Related publications 



VvWRKY2 [29,30] 



VvWRKYl [18] 



VpWRKY2 [30] 



VpWRKYl [30] 
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Table 1 Identified WRKY genes in 12x V. vinifera 'Pinot Noir' genome (Continued) 



GSVI\^01 032661 001 


VvWRKY45 


III 


ehrl3 


289 




GSVI\^0 1032662001 


VvWRKY46 


III 


ehrl3 


309 




GSVI\^0 1033063001 


VvWRKY47 


lie 


ehrl4 


183 




GSVI\^01 0331 88001 


VvWRKY48 


lid 


chr4 


268 


WRKYl 1 [33,34] 


GSVI\^01 0331 94001 


VvWRKY49 


lie 


chr4 


157 




GSVI\^01 0331 95001 


VvWRKYSO 


lie 


chr4 


102 




GSVI\^01 0341 48001 


VvWRKYSl 


lie 


chrB 


300 




GSVI\^0 1034968001 


VvWRKY52 


lie 


chr5 


310 




GSVI\^0 1035426001 


VvWRKY53 


lie 


chr4 


167 




GSVI\^01 035884001 


VvWRKY54 


lla 


chr4 


263 




GSVI\^01 035885001 


VvWRKY55 


lla 


chr4 


287 


VpWRKYS [32] 


GSVI\^01 035965001 


VvWRKY56 


1 


chr4 


531 




GSVI\^0 1036223001 


VvWRKY57 


lid 


ehrl4 


305 




GSVI\^0 1037686001 


VvWRKY58 


lib 


ehrl9 


497 




GSVI\^01 037775001 


VvWRKY59 


1 


ehrl9 


553 





Categorization of VvWRKYs basis on conserved WRKY 
domains 

The disposition of structural domains in amino acid se- 
quences is an important clue to analyze the evolution 
and relationship between highly divergent sequences 
[38]. The relationships among the 59 WRKY proteins 
were investigated through constructing phylogenetic 
trees based on multiple alignments of the predicted 
amino acid sequences of the WRKY domains. As shown 
in Figure 1, we classified the 59 VvWRKY proteins into 
four large groups according to the results of the phylo- 
genetic analyses. The models of conserved amino acid 
sequences of WRKY domain and zinc-finger structure in 
four groups were shown in Additional file 2: Figure S2. 
Twelve of the WRKY proteins contained two complete 
WRKY domains and a C2H2-type zinc finger motif. 
These proteins constituted group I. The N-terminal 
WRKY domain (NTWD) and C-terminal WRKY domain 
(CTWD) of VvWRKY27, VvWRKY41 and VvWRKY56 
were clustered into a same clade in group I. According 
to Eulgem et al. [9] and by using WRKY proteins in Ara- 
bidopsis as references, 39 VvWRKY in group II were cat- 
egorized into five subgroups. Three members were 
found in subgroup Ila, 8 in lib, 16 in lie, 6 in Ild and 6 
in He. Group II was divided into two parts. Subgroup 
Ila, lib and lie showed a close relationship with Group 
III WRKY proteins. And subgroups lid and He belonged 
to a separate clade which was closely related to group 
IV. Subgroup He showed higher divergence than the 
other subgroups. There were also 6 WRKY proteins in 
group HI, and 2 in group IV which lacked a complete 
zinc-finger structure. 



RT-PCR based transcription levels detection of VvWRKYs 
in different tissues 

To investigate if the putative VvWRKYs were expressed 
and assess their transcription levels in grape, we exam- 
ined the expression of these genes in different grape tis- 
sues. Among all VvWRKYs ^ we successfully designed and 
verified 58 primer pairs representing all candidate 
VvWRKYs except for VvWRKYSS (Figure 2). All tran- 
scripts can be detected at least in one tissue. Nineteen 
VvWRKYs (including VvWRKY02, 11, 12, 13, 14, 17, 20, 
24, 28, 30, 33, 34, 3S, 36, 39, 41, 42, 48 and 52) were 
found expressed in all tissues used. Six VvWRKYs 
{VvWRKYOS, 09, 22, 40, 44 and 58) were found only 
expressed in young tissues. VvWRKYOS was expressed in 
the stem apex and young fruit. VvWRKY40 was found in 
stem apex, young fruit and root. VvWRKY09, 22, 44 and 
58 were detected in young leaf, stem apex, young fruit 
and root. 

Gene-chip based expression analysis of 26 VvWRKYs 
under various stresses 

Although we identified WRKY transcription factors from 
the V vinifera genome, functions for these genes in re- 
sponse to abiotic and biotic stress remain unknown. 
Using microarray results from publically available data, 
it was possible to find gene expression data from mul- 
tiple experimental conditions for several of the grapevine 
WRKY genes. We carefully checked the genes on the 
'GeneChip Vitis vinifera (Grape) Genome Array (Affy- 
metrix) and a total of 26 VvWRKYs were found on this 
chip. Microarray data related to salinity, water-deficit, 
PEG, cold, ABA and pathogen stresses were downloaded 
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Figure 1 Phylogenetic tree of VvWRKYs. The unrooted phylogenetic tree of WRKY domains was constructed with MEGA5.1 program with the 

neighbor-joining method. The numbers beside the branches represent bootstrap values based on 1000 replications. The name of groups (I, II, III and 

IV) and subgroup (a-e) were shown at the outside of the circle. The WRKY named with suffix -N or -C indicated the N-terminal WRKY domain (NTWD) 

or the C-terminal WRKY domain (CTWD) in one VvWRKY with two WRKY domains. AtWRKYs were used as reference to categorize VvWRKYs. 
J 



and their corresponding probes and the CV (coefficent 
of variation of the corresponding treatment means) of 
these genes in each of the microarray experiments were 
listed in Additional file 3: Table SI. If the expression of a 
probe set (gene) is affected by some of the treatments in 
an experiment, it shows a higher CV (more fluctuation); 
and vice versa. According to the data, the CV of 20 of 
the 26 VvWRKYs were over 5% in at least one experi- 
ment. The highest CV appeared in VvWRKY57 (up to 
36%) associated with compatible viral diseases in berry 
experiment in V. vinifera cv. Cabernet Sauvignon. 
VvWRKYOS, 06, 08, 28, and 55 responded to both abiotic 
and pathogens stresses while VvWRKY21, 39, 48 seemed 
to respond primarily to pathogens stresses. 

To test the correlation between the expression pat- 
terns of 26 VvWRKYs and their phylogenetic relation- 
ship, a hierarchical cluster analysis was performed using 
the 11 stress related experimental datasets (Figure 3). 
Red, black and green elements in the matrix indicate up-, 
no change- and down-regulated expression of WRKY 
transcription factors, respectively. From the heat map, 
twenty-six genes were clustered into four clades. Care- 
fully analyzing the cluster of expression data in response 
to abiotic stresses experiments and comparing this with 
the VvWRKYs phylogenetic tree, we found that genes 
with close phylogenetic relationship were classified into 



the same clade during hierarchical cluster analysis. The 
most obvious evidence can be found in clades 3 with 5 
WRKY subgroup Ild genes (including VvWRKY07, 24, 
39, 48 and 57), which show similar expression patterns 
in response to salt, PEG and cold stresses. Clade 1 con- 
tained three WRKY group I genes and two group IIC 
genes. Clade 2 was mainly composed by WRKY group I 
and IIC and contains a majority of cold stress-related 
VvWRKYs (Also shown in Additional file 4: Table S2). 
Clade 4 only had one gene and that gene was from 
WRKY group III. 

Real-time RT-PCR based expression analysis of VvWRKYs 
under cold treatment in V. vinifera 

To examine the response of VvWRKYs under cold stress 
in grape, we examined the transcription levels of 
VvWRKYs in shoot apices of 'Muscat Hamburg' under 
cold-treatment (4°C). VvWRKYOS, 21, 32 and 40 were 
excluded from cold-treated experiment since their Ct 
value of amplification curve were over 35 cycles in the 
templates of normal and cold-treated shoot apex. De- 
tected VvWRKYs can be classified into four groups ac- 
cording to expression patterns as shown in Figure 4 and 
Additional file 5: Figure S3: A) sustained up-regulated 
during cold treatment (22 genes. Figure 4A), B) changed 
above 2 fold with irregular pattern (9 genes. Figure 4B), 
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YL MLT S YF RF I 



YL MLT S YF RF R 



VvWRKYOl 
VvWRKYOZ 
VvWRKYOS 
VVWRKY04 
VvWRKYOS 
VvWRKYOG 
VVWRKY07 
VvWRKYOS 
VVWRKY09 
VvWRKYlO 
VvWRKYll 
VvWRKYlZ 
VvWRKYlS 
VVWRKY14 
VvWRKYlS 
VvWRKYie 
VvWRKYl? 
VvWRKYlS 
VVWRKY19 
VvWRKYZO 
VvWRKYll 
VVWRKY22 
VVWRKY23 
VVWRKY24 
VVWRKY25 
VVWRKY26 
VVWRKY27 
VVWRKY2S 
VVWRKY29 
VvMDH 



VvWRKYSO 
VvWRKYSl 
VVWRKY32 
VVWRKY33 
VVWRKY34 
VVWRKY35 
VVWRKY36 
VVWRKY37 
VVWRKY39 
VVWRKY40 
VvWRKY41 
VVWRKY42 
VVWRKY43 
VVWRKY44 
VVWRKY45 
VVWRKY46 
VVWRKY47 
VVWRKY4S 
VVWRKY49 
VvWRKYSO 
VvWRKYSl 
VVWRKY52 
VVWRKY53 
VVWRKY54 
VVWRKY55 
VVWRKY56 
VVWRKY57 
VvWRKYSS 
VVWRKY59 
VvACT 



Figure 2 RT-PCR analyses of presence of WI/l//?/C/ transcripts in seven grape tissues. YL: young leaf; ML: mature leaf; T: tendril; S: stem apex; 
YF: young fruit; RF: ripe fruit; R: root. VvMDH and VvAG' were used as control. 



C) sustained down-regulated (5 genes, Figure 4C) and D) 
no significant difference (18 genes, as shown in Additional 
file 5: Figure S3). The relative expressions of 36 genes 
(Figure 4A, B and C) were significantly different as cold 
treatment. The greatest increase in expression (nearly 30 
fold) was found in VvWRKYSS at 48 h cold treatment. 
VvWRKYlS and VvWRKY46 had the largest up-regulation 
of greater than 6 fold at 8 hours after cold treatment. 
While VvWRKYlS was degraded after 24 hours, the ex- 
pression of VvWRKY46 demonstrated both up and down 
regulated with a spike of expression at 48 hours after in- 
tensive degradation at 24 hours. 

Exogenous ABA induced accumulation of VvWRKYs in 
V, Vinifera 

To illustrate how the VvWRKYs respond to ABA and 
whether the cold stress related VvWRKYs may partici- 
pate in the ABA-dependent cold signal pathway, ABA 



treated grapevine apices were examined using qRT-PCR. 
VvWRKY12, 29, and 46 were excluded from this ex- 
periment due to their higher Ct value (Figure 4D and 
Additional file 6: Figure S4). Among the 55 VvWRKYs 
we detected, twelve VvWRKYs were expressed over 2- 
fold greater within 2 h of exogenous ABA treatment 
(Figure 4D). After statistical analyses of qRT-PCR re- 
sults, 7 of them were evaluated to significantly change 
during exogenous ABA treatments. Transcripts of 
VvWRKYSS showed the greatest increase in expression 
at 0.5 h after ABA treatments. Six other genes showed 
increases in expression 1 h after exogenous ABA treat- 
ment (Figure 4D). 

When the data from the cold and ABA experiments 
were compared, 6 of 7 genes {VvWRKY, 19, 2S, 3S, 42, 
SO and SS) that were up-regulated during exogenous 
ABA treatment were also up-regulated under cold treat- 
ment (Figure 4A and B, marked by underline). Two 
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Figure 3 Cluster analyses of VvWRKYs from 16 k Affymetrix V. 
vinifera gene-chip data in PLEXdb database. The relative 
expression values of 26 VvWRKYs responding to different abiotic 
stresses (salinity, water deficit, cold) were used in analysis. Red, black 
and green elements in the matrix indicate up-regulated, no change 
and down-regulated WRKY genes, respectively. Those genes can be 
classified into four groups according to expression patterns, which 
were shown in different color with its group IDs that coincide with 
Figure 1. The red color was used to emphasize the VvWRKYs that 
changed in expression over 2 fold under cold treatment. 



genes {VvWRKYSS, 28) were greatly up-regulated, over 
10 fold. The expression levels of the rest of the 44 
VvWRKYs were lower than 2-fold and not significantly 
changed during exogenous ABA treatments (Additional 
file 6: Figure S4). 

Identification of candidate cold-stress related VvWRKYs 

Previously we reported the changes of the transcriptome 
during cold-treatments in 'Muscat Hamburg' and identi- 
fied 14 cold-stress related VvWRKYs (we reported 16 
VvWRKYs but subsequent annotation of these genes 
allowed us to exclude two genes that do not belong to 
the WRKY gene family) [36]. Gene-chip based methods 
also allowed to identify 10 cold-stress related VvWRKYs 
[39]. In order to overcome the deficiencies of determin- 
ing gene expression from a single technological approach 
and obtain more reliable results, we compared the data 
from three different methods. Fourteen VvWRKYs from 
our previous transcriptome analysis, ten from publically 
available gene-chip based data and 36 genes from qRT- 
PCR results (this study) were used. The results were sum- 
marized in Figure 5 and Additional file 4: Table S2. Three 
VvWRKYs {VvWRKY12, 28, 55) showed identical expres- 
sion patterns and were found up-regulated over 10 fold in 
at least one time-point under cold-treatment by qRT-PCR 



(Figure 4A). A total of 12 VvWRKYs were confirmed by 
two experimental methods (Figure 5 A and B). VvWRKYS6 
was identified as up-regulated gene under cold treatment 
only in the gene-chip studies. Twenty-two genes that were 
characterized by qRT-PCR were not supported by the 
other studies. It is worth mentioning that down-regulated 
VvWRKYs under cold-treatment were only identified by 
qRT-PCR based method. 

Discussion 

WRKY family in grape 

Considering the important roles that WRKY transcrip- 
tion factors play during plant development and in re- 
sponse to various stresses, it is not surprising that we 
identified so many family members in grapevine. Previ- 
ously, 74 WRKY genes were found in Arabidopsis [2], 55 
in cucumber [40], 102 in rice [2], 47 in castor bean [41], 
86 in Brachypodium distachyon [42] and 136 in maize 
[43]. Here we identified 59 candidate WRKY proteins in 
V, vinifera and categorized them into four groups. 

Group I WRKY proteins 

When compared with WRKY family groups, WRKYs in 
primitive plant ancestors Giatdia lamhlia, Dictyostelium 
discoideum and Chlamydomonas reinhardtii closely re- 
sembled Vitis group I [7,38]. In our study, two domains 
of VvWRKYs in group 1 were closely related. A BLASTp 
search of EuGene.l 100010359 from an ancient alga spe- 
cies {Ostreococcus sp, RCC809) which has a single WRKY 
domain allowed us to identify 9 corresponding WRKY ho- 
mologs in grape and 8 of these belonged to group I by 
MAP VIEW (Plant Genome Duplication Database) [44]. 
These data support the hypothesis that the dual WRKY 
domains present in members of group I may be derived 
from a single WRKY domain duplication [6,7]. 

Group II WRKY proteins 

Group II was divided into three parts: subgroup Ila + lib, 
subgroup lie and subgroup Ild + He (Figure 1). Subgroup 
Ila + lib belong to the same clade and is sister to the 
WRKYs in group I. Interestingly, the presumed function of 
CTWDs in group I for sequence- specific DNA binding [9] 
were more similar to the single WRKY domain members 
in group II and III than to the NTWDs of group I, This re- 
sult may indicate that subgroup Ila + lib evolved from 
group I by domain structure loss of the group I NTWD. 

Group III WRKY proteins 

Group III in the phylogenetic tree was most closely related 
to the very large subgroup lie, which was separately into 
four clades and seemed to indicate an expansion of the 
gene family. A thorough search of the Plant Transcription 
Factor Database (http://planttfdb.cbi.pku.edu.cn) indicated 
that the earliest evolutionary occurrences of group III 
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(See figure on previous page.) 

Figure 4 qRT-PCR assays of the expression patterns of VvWRKYs under cold and exogenous ABA treatments. Tlie default expression 
value for each gene was 1 at 0 hours before treatment. A, B and C represent the subgroups with different expression patterns in cold treatment 
and D represents the genes that up-regulated over 2 fold in ABA treatment. A: sustained up-regulated genes in cold treatment; B: genes that 
changes over 2 fold but without significant tendency in cold treatment; C: sustained down-regulated genes in cold treatment; D: up-regulated 
genes in exogenous ABA treatment. VvWRKYs that accumulated in both cold and exogenous ABA treatments were underlined. One-Way ANOVA 
analysis was used to test the impact of timing of cold treatment. When the effects were significantly different, we examined the difference 
between treatments using post hoc multiple comparisons (LSD, p < 0.05). All data analyses were conducted using IBM SPSS Statistics 20, and the 
results were displayed through a, b, c and d. 



genes were those found in ferns {Selaginella moellendorf- 
fii). There was no evidence of any sequenced plant species 
that only contain members of group I and III but we 
found in some species with only members of group I and 
II, for example in mosses {Physcomitrella patens) [1], and 
some gymnosperms (Pinus taeda) (http://planttfdb.cbi. 
pku.edu.cn). We speculated that group III may have 
evolved from group II, particularly lie. As group III 
WRKYs in Arabidopsis responded to diverse bio tic stresses 
[45], group III members may indicate adaptation of early 
plants to the stressful conditons associated with the 
colonization of land and subsequent increase in biotic 
pathogen pressures. 

Group IV WRKY proteins 

We found that group IV WRKY proteins, which were 
characterized by the loss of the zinc-finger domain, were 
in the same clade as subgroups Ild + He. VvWRKY02 
and Vv WRKYS 7 were duplicated gene pairs according to 
a whole genome analysis of grapevine gene duplications 
[46]. This might suggest an origin of group IV from sub- 
groups lid + He. Group IV proteins were considered 
non-functional due to the loss of the zinc-finger domain 
[10]. However, these genes of group IV can be found in 
all higher plant species as well as in algaes (Bathycoccus 



prasinos: Bathyl7g02050). Furthermore, some genes 
were expressed in rice {OsWRKYS6) [10] as were two genes 
identified in this study {VvWRKY02 and VvWRKY29). 
Therefore, it remains questionable whether group IV 
WRKYs have biological function in plants. 

VvWRKYs participate in development and stress-related 
signal pathways 

WRKYgQViQs were found to be expressed in many tissues 
and seem to be involved in regulating plant developmen- 
tal and physiological processes. Transcriptomic analysis 
of senescence in the flag leaf of wheat demonstrated that 
W7?/<rF transcription factors are greatly up-regulated dur- 
ing the senescence process [47]. OsWRKY78 was found 
to be up-regulated in elongating stems and knockdown 
mutations in this gene cause plants to produce a semi- 
dwarf and small seed phenotype caused by reducing cell 
length [48]. Moreover, the transcription of GhWRKYlS 
was observed abundant in the roots and stems of to- 
bacco and transgenic overexpression lines of these plants 
displayed faster elongation at the earlier shooting stages 
[49]. Here the expression of 15 VvWRKYs (Figure 2) can be 
detected in all grape tissues we used, which may indicate 
its fundamental roles in different cell-types in grape. Simi- 
lar to expression patterns observed in other plant species. 




Figure 5 An overview of cold stress-related VvWRKYs in three sets of data. A: The Venn diagram of the cold stress-related VvWRKYs obtained 
from qRT-PCR, transcriptome and gene-chip data. B: The VvWRKYs that were found in more than one type of experimental data. Green color in 
forms indicated VvWRKYs induced by exogenous ABA. 

V ) 



Wang ef al. BMC Plant Biology 2014, 14:103 
http://www.bionnedcentral.conn/1471-2229/14/103 



Page 10 of 14 



VvWRKYs were found to be expressed in young tissues 
such as young leaf, shoot apex, tendril and young fruit. 

Several numbers of VvWRKYs were found activated in 
more than one type of stress condition (Figure 3 and 
also Additional file 3: Table SI). VpWRI<Y3, homologous 
to VvWRKYSS was observed to be up-regulated in re- 
sponse to many different sources of stress, including 
pathogen exposure, salicylic acid, ethylene, cold and 
drought stress [32]. VvWRKYs that were up-regulated in 
response to more than two types of stresses (e.g. patho- 
gen and drought) supported the occurrence of cross-talk 
between signal transduction pathways in response to dif- 
ferent stress conditions in plants [50]. 

Phylogenetic relationships between VvWRKY genes 
suggested that there may be conserved responses of these 
genes to salt exposure, PEG and cold-stress (Figure 3). All 
members of group Ild clustered into one clade with simi- 
lar expression pattern during these three stress conditions, 
suggesting the function of these VvWRKY proteins may 
relate to the structures of WRKY domains. Subgroup lid 
was identified as a novel CaM-binding transcription factor 
family in plants and their conserved structural motif was a 
Ca^"^-dependent CaM-binding domain [51]. Thus the 
placement of the WRKYs in the phylogenetic tree may 
also help to predict function of new members that belong 
to certain gene family. 

VvWRKYs that participate in the cold related signal 
transduction in grape 

Three different experimental methods were combined to 
robustly analyze the response of VvWRKY genes to cold 
stress (Figure 5 and Additional file 4: Table S2). Results 
from qRT-PCR demonstrated the greatest number of 
cold stress-related VvWRKYs (36) while gene-chip based 
methods identified the least, 10 VvWRKYs, This differ- 
ence may be attributed to the method used but is also 
likely due to differences in the treatment conditions be- 
tween experiments. During Digital Gene Expression pro- 
file (DGE) analysis [36], plant material was obtained 
from 4 h cold treatment at 4°C, whereas in our pRT- 
PCR experiment, we used samples collected at several 
different time periods (at 8 h, 24 h and 48 h after cold 
treatment at 4°C). Additionally, multiple matched tags 
were excluded from the final analysis performed by Xin 
et al. [36], which may have reduced the number of iden- 
tified cold related VvWRKYs. Finally, gene-chip based 
methods may bias results due to a lower number of 
genes with corresponding probes related to the WRKY 
proteins (only 26 WRKY). By integrating the data from 
different methods, we obtained more reliable results and a 
total of 15 candidate cold tolerance VvWRKYs (Figure 5) 
were identified during our investigation. 

According to previous studies, the transcriptional con- 
trol of plant responses to cold stress can be divided into 



ABA-dependent and ABA- independent signal pathways 
[52]. The results of our study also indicated that 15 pu- 
tative cold stress-related VvWRKYs can be divided into 
two groups according to their responses to exogenous 
ABA. Three VvWRKYs {VvWRKY28, 42 and 55) may par- 
ticipate in an ABA-dependent signal pathway and other 
12 in ABA-independent pathway. WRKY transcription 
factors have been identified as key components in the 
ABA signaling pathways [8,53]. In rice, OsWRKY24, 51, 71 
and 72 are induced by (ABA) in aleurone cells. 
OsWRKY24 and 45 were functional as negative regula- 
tors in ABA induction of the HVA22 promoter-beta- 
glucuronidase construct, while OsWRKY72 and 77 
synergistically interacted with ABA to activate this 
reporter construct [10]. It is still unknown how 
WRKYs participate in the cold stress -related signal 
pathway and what relationship these genes have with 
C-repeat Binding Factor genes (CBFs), which are critical 
transcription factors responsible for cold tolerance in 
plant [54]. 

The reliability of the identified 15 cold-related 
VvWRKYs was also supported by homologous genes in 
other species. STHP-64, which showed high similarity 
with VvWRKY43, was not present in leaves until Novem- 
ber and December in Solanum dulcamara [55]. WRKY38, 
a homolog gene of VvWRKY14, was transiently accumu- 
lated when leaves and roots were exposure to low 
temperature in barley [56]. BcWRKY46 showed higher 
similarity with VvWRKYSS and responded to low temper- 
atures in Pak-choi. Constitutive expression of BcWRKY46 
reduced the freezing susceptibility in transgenic tobacco 
[57]. The transcription level of VvWRKY55 was up- 
regulated robust under cold treatment. Its homolog gene, 
WRKY71 was found in banana with a similar expression 
pattern [58]. All these VvWRKYs mentioned above were 
confirmed by at least two set of experiment methods, 
which provided appropriate candidates to illustrate the 
roles of WRKY protein under low temperature-related sig- 
nal pathways in grape. 

Although low-temperature related WRKYs were iso- 
lated in several species, the mechanism of how WRKYs 
respond to cold signals and regulate the expression of 
downstream genes is still largely unknown. Further work 
is needed to elucidate the function of these important 
genes in low-temperature related signal pathways. Previ- 
ously we reported the different expression patterns of 
WRKYs in V amurensis, a cold-hardness species. The 
WRKY genes identified here from V vinifera may accel- 
erate the functional analysis of this gene family in V, 
amurensis. The comprehensive analysis of cold stress- 
related WRKYs in two different Vitis species with con- 
trasting cold hardiness phenotypes would certainly help 
to illustrate the function of WRKY genes in conveying 
cold hardiness in grapevine. 
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Conclusions 

In summary, a total of 59 VvWRKYs in the V. vinifera 
genome were identified. The VvWRKYs were unevenly 
distributed in 18 of the 19 chromosomes. WRKY do- 
main based phylogenetic analysis allowed categorizing 
59 VvWRKYs into four large groups. A majority of 
VvWRKYs were found expressed in more than one tissue 
in V, vinifera. Gene-chip based data analysis suggested 
that a subset of VvWRKYs was activated in respond to 
diverse biotic and abiotic stresses. The transcription level 
of 36 VvWRKY genes changed over 2 fold after cold in- 
duction. A comparative analysis of qRT-PCR results, 
gene-chip based data and transcriptome analysis allowed 
us to identify 15 VvWRKYs that show identical expres- 
sion patterns during cold treatment at least in two kinds 
of analyses. These studies not only increase our know- 
ledge of WRKY family, but also provide candidate genes 
for future functional analysis of VvWRKYs involved in 
the low temperature-related signal pathways in grape. 

Methods 

Identification of WRKY genes in the grape genome 

Candidate WRKY proteins were identified from the 12X 
V, vinifera cv. Pinot noir genome (quasi-homozygous 
line PN40024, http://www.phytozome.net). Full-length 
amino acids sequences of all WRKY proteins in Arahi- 
dopsis thaliana (http://www.arabidopsis.org/) were used 
as query sequences. A BLASTp search was performed 
and £- value of e"^ was used as the threshold [59]. Candi- 
date WRKY proteins were manually confirmed [60] by 
searching for WRKY domains in the candidate amino 
acids sequences using SWISS-MODEL (http://swissmo- 
del.expasy.org/) and the results were shown in Table 1. 

Phylogenetic analysis of WRKY family 

Multiple alignments of the amino acid sequences of 73 
WRKY domains from V, vinifera were performed using 
CLUSTALW by MEGA5.1 [61]. Twelve Arabidopsis 
WRKY domains from different WRKY groups were used 
as references to categorize the WRKY proteins from 
grape. The GenBank accession numbers of those AtWR- 
KYs are AtWRKYOl: ABJ17102, AtWRKYll: AEE85928.1, 
AtWRKY14: AAP21276.1, AtWRKYlS: AAM78067, 
AtWRKY21: AAB63078.1, AtWRKY27: ABH04558, 
AtWRKY28: AEE84006, AtWRKYSl: AEE84546.1, 
AtWRKY38: AED93044.1, AtWRKY41: AEE82969, 
AtWRKY43: AEC10646.1, AtWRKY45: ABD57509.1, 
AtWRKY49: AAQ62425.1. The parameters used during 
alignment were: protein weight matrix: Gonnet series; 
negative matrix: on; gap open penalty: 10; gap extension 
penalty: 0.20; delay divergent sequences: 30; residue- 
specific gap penalties: on; hydrophilic penalties: on; gap 
separation distance: 0; end gap separation penalty: on. 
An unrooted phylogenetic tree was constructed using 



Neighbor-Joining (NJ) methods and bootstrapped with 
1,000 iterations to help identify WRKY protein groups. 

Plant materials 

'Muscat Hamburg' {V, vinifera) was obtained from the 
Institute of Botany, the Chinese Academy of Sciences. 
Tissues of young leaves, mature leaves, tendril, stem 
apex, root, young fruits and ripe fruits were collected 
from the vineyard in July, 2012. Cold and exogenous 
ABA treatment experiments were performed on tissue 
culture seedlings of 'Muscat Hamburg' according to Li 
et al. [62]. Briefly, seedlings were cultured on 1/2 B5 
medium with 30 g/L of sucrose in a growth chamber 
under 16-h light/8-h dark photoperiod at 26°C. Cold 
treatments were performed in another growth chamber 
with the same parameters except for temperature (4°C). 
Seedlings with five well developed leaves were used and 
the shoot apex with one well developed leaf was col- 
lected at 0 hour (h, used as control), 8 h, 24 h and 48 h. 
Seedlings with five well developed leaves were trans- 
planted in 1/2 B5 nutrient solution. Exogenous ABA 
treatments were performed after one week under normal 
culture conditions with 100 (iM ABA and the shoot apex 
with one well developed leaves were collected at 0 h 
(used as control), 0.5 h, 1 h and 2 h after treatments. 
Three independent replicates were collected for each 
time point and frozen in liquid nitrogen. Samples were 
then stored at - 80°C for the following RNA isolation. 

Expression patterns analysis of VvWRKYs by RT-PCR 

Total RNA was isolated from collected samples using 
Plant Total RNA Isolation kit (Tiandz Inc; Beijing, 
China). RNase-free DNase (RQl, Promega) was used to 
degrade DNA from total RNA. cDNA was synthesized 
by the Superscript III Reverse Transcriptase (Invitrogen) 
with 01igo(dT)i8 (Promega) according to the manu- 
facturers instructions. Primer pairs (Additional file 7: 
Table S3) for VvWRKYs were designed by Primer 3 
(http://bioinfo.ut.ee/primer3-0.4.0/) and tested by NCBI 
Primer BLAST. Two genes, p-actin (GenBank accession: 
EC969944; sense primer: 5'-CTTGCATCCCTCAGC 
ACCTT-3'; antisense primer: 5'-TCCTGTGGACAAT 
GGATGGA-3') and malate dehydrogenase gene (MDH; 
GenBank accession: EC921711; sense primer: 5'-CCAT 
GCATCACCCACAA-3'; antisense primer: 5'-GTCAA 
CCATGCTACTGTCAAAACC-3') were used as positive 
control for RT-CR [63]. Three biological replicate and 
35 cycles for each reaction were performed. PCR prod- 
ucts were detected by agarose gel electrophoresis with 
2.5% gel concentration. 

Gene-chip based expression pattern analysis of VvWRKYs 

We explored the expression profiles of VvWRKYs using 
publically available data from the 16 k Affymetrix V, 
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vinifera gene-chip stored at PLEXdb (Plant Expression 
Database) [64] to explore the response of VvWRKYs dur- 
ing biotic and abiotic stresses in grape. The different 
studies and datasets that were included in these analyses 
were: A) a short term abiotic stress experiment in 'Cab- 
ernet Sauvignon' [39], B) a long-term salt and water 
stress study [65]; C) a study examining gene expression 
associated with compatible viral diseases in grapevine 
cultivars [66]; D) an experiment designed to examine the 
powdery mildew-induced transcriptome in a susceptible 
grapevine 'Cabernet Sauvignon' [67]; E) the com- 
plimentary dataset of the powdery mildew-induced tran- 
scriptome of a resistant grapevine 'Norton' [67]; F) a 
study of gene expression in grapevine in response to Bois 
noir infection [68]; G) a study of the grape skin tran- 
scriptome of berries grown on an exogenous abscisic 
acid treated vine [69]; H) the complimentary dataset of 
the grape skin transcriptome in the berries cultured 
in vitro and treated with exogenous ABA [69]; and lastly, 
I) a gene expression study associated with compatible 
viral diseases in the berry [70]. In our comparative ana- 
lysis, we divided these experiments into either abiotic or 
biotic stresses related datasets. For each microarray ex- 
periment, the Affymetrix MAS5.0 normalized data were 
used for calculations of the fold change of differentially 
expressed genes. Probe sets corresponding to the puta- 
tive VvWRKYs were identified at PLANEX (http://pla- 
nex.plantbioinformatics.org) and completed via PLEXdb 
blast tool. Comparisons of WRKY expression level from 
gene-chip data for the short term abiotic stress treat- 
ment in 'Cabernet Sauvignon' was performed using 
Cluster 3.0 and JavaTreeview. 

Quantitative RT-PCR 

Total RNA was isolated from cold and exogenous ABA 
treated shoot apices following the cDNA synthesis 
methods mentioned above. Synthesized cDNA was di- 
luted 1:10 with ddH20, and the quantitative RT-PCR re- 
action mixture contained 5 (il of 2 x SYBR Green I 
Master Mix (Roche, USA), 2.6 \xL ddH20, 0.2 \xL of 
10 (iM solution of each primer and 2 (iL diluted tem- 
plate cDNA. Reaction specificities for each primer pair 
was tested using qRT-PCR melting curve analysis. The 
experiment was carried out using a StepOnePlus real- 
time PCR Instrument (Applied Biosystems). Transcrip- 
tion levels of each VvWRKY was normalized against 
the average of |3-actin, MDH (as mentioned above) 
and glyceraldehyde-3-phosphatedehydrogenase (GAPDH: 
CB973647; sense primer: 5'-TTCTCGTTGAGGGCT 
ATTCCA-3'; antisense primer: 5'-CCACAGACTTCAT 
CGGTGACA-3') [63]. Each sample had three biological 
and two technical replicates to ensure the accuracy of re- 
sults, and RNA samples with the same reverse- 
transcription (without Reverse Transcriptase) and dilution 



procedure were used as negative controls. The Ct values 
and the real-time PCR efficiencies were obtained using 
Lin-RegPCR [71] and the normalized relative quantities 
and standard errors for each sample were calculated by 
qbaseplus [72]. The relative expression level of each 
VvWRKY in different templates was calculated based on 
normalized relative quantities. We used One- Way ANOVA 
analysis to test the impact of timing of cold treatment. 
When the effects were significantly different, we examined 
the difference between treatments using post hoc multiple 
comparisons (LSD, p < 0.05). All data analyses were con- 
ducted using IBM SPSS Statistics 20. 

Additional files 



Additional file 1: Figure SI. Chromosomal location of 57 VvWRKYs. 
VvWRKY03 was located on 'chromosome 1 random' and VvWRKY04 was 
located on 'chromosome unknown'. Neither was shown here. 

Additional file 2: Figure S2. The models of conserved amino acid 
sequences of WRKY domain and zinc-finger structure in four groups. The 
numbers behind the charts indicated gene numbers in each group. 

Additional file 3: Table SI. The coefficient of variation of the 
corresponding treatment means (CV) and probe set IDs of VvWRKYs in 9 
experiments. A higher CV means the expression of the probe set (gene) 
is affected by treatments in an experiment. Five VvWRKYs that didn't 
show any changes in any treatments are marked by green color. 

Additional file 4: Table S2. Cold stress-related VvWRKYs obtained in 
one of three experimental methods. Yellow, red and blue forms represent 
genes obtained via qRT-PCR, gene-chip data and transcriptome data 
respectively. Exogenous ABA induced VvWRKYs were shown in green 
color in form. 

Additional file 5: Figure S3. Quantitative RT-PCR assays of the expression 
level of 18 VvWRKYs under cold treatment. The transcription level of these 
genes didn't show significant changes during cold treatment in V. vinifera. 

Additional file 6: Figure S4. Quantitative RT-PCR assays of the expression 
patterns of 44 VvWRKYs under exogenous ABA treatment. The transcription 
level of these genes didn't show significant changes during exogenous ABA 
treatment in V. vinifera. 

Additional file 7: Table S3. The primers used for expression pattern 
analysis for VvWRKYs. 
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